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Circular microresonators are micron sized dielectric disks embedded in material of lower refractive 
index. They possess modes of extremely high Q-factors (low lasing thresholds) which makes them 
ideal candidates for the realization of miniature laser sources. They have, however, the disadvantage 
of isotropic light emission caused by the rotational symmetry of the system. In order to obtain high 
directivity of the emission while retaining high Q-factors, we consider a microdisk with a pointlike 
scatterer placed off-center inside of the disk. We calculate the resulting resonant modes and show 
that some of them possess both of the desired characteristics. The emission is predominantly in the 
direction opposite to the scatterer. We show that classical ray optics is a useful guide to optimizing 
the design parameters of this system. We further find that exceptional points in the resonance 
spectrum influence how complex resonance wavenumbers change if system parameters are varied. 

PACS numbers: 42.55.Sa, 42.25.-p, 42.60.Da, 05.45.Mt 



I. INTRODUCTION 

Thin dielectric microcavities of different shapes are 
widely used as microresonators in laser physics and in- 
tegrated optics (see [TJ [2] and references therein). The 
new directions in microcavity laser research have been re- 
cently reviewed in [3^ . Microresonators are open systems 
coupled to the external world, and therefore do not have 
bound states. Instead, their eigenmodes (resonances) are 
characterized by complex wavenumbers k = kr + iki = 
ijjjc — i/{cT). Here r is the lifetime of the resonance 
and c denotes the speed of light. The so-called reso- 
nance quality factor {Q factor for short) is then defined 
as Q — kr/2\ki\, and is a measure of the suitability of a 
mode for lasing. 

The resonances of any passive (no active lasing parti- 
cles) microcavity filled with nonmagnetic (B = H) dielec- 
tric material are the solutions of the time-independent 
Maxwell equations 



VxE^ifcH, VxH = -ifce(r)E 



(1) 



when appropriate boundary conditions on the fields H 
and E are imposed. When the fields are independent of 
z, which is strictly speaking only the case for a cylindri- 
cal microcavity of infinite length (see below) , the above 
equations in cylindrical coordinates (r, (p, z) reduce to 
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The solutions can be separated into TM waves 'trans- 
verse magnetic, i.e. — 0) and TE waves 'transverse 
electric', i.e. E^ — 0). Introducing the position depen- 
dent refractive index n{r, ip) = y^£(r, p), we get for TM 
modes a scalar wave equation for Ez, 
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The magnetic field El is then recovered from 
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For TE modes, it is advantageous to introduce the new 
field hz{r, p) = i?z(r, p) /n{r, p). This field satisfies again 
a scalar wave equaton. 
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and the electric field E can be recovered in this case from 



Er — 



krn{r, p) 
i 



kn{r,p) 



dhz 1 dn 

dp n{r, p) dp 

dhz 1 dn 

— - H 

dr n{r, p) dr 



hz{r, p) 
hz{r, p) 



(7) 



The applicability of these equations can be extended from 
infinitely long cylindrical microcavities to flat microcav- 
ities for which the cavity thickness is only a fraction of 
the mode wavelength. The refractive index has then to be 
replaced by an effective refractive index rioff(r, 93) which 
takes into account the material as well as the thickness 
of the cavity. We note that it is in general a difficult 
problem to explicitly compute nf.f{{r,p) for a given ma- 
terial and thickness (see, e.g., Ref. U), and one therefore 
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typically resorts to using an experimentally determined 
effective refractive index instead. 

The simplest microcavity shape is a thin circular mi- 
crodisk. The technological progress in recent years has 
made possible the construction of microdisks in the /im- 
domain. They are natural candidates for the construction 
of lasers since some of their modes have extremely high 
Q- factor (low thresholds) [5] [5]. In those modes, which 
are called 'whispering gallery modes', light is trapped by 
total internal reflection and circulates along the circum- 
ference of the disk. As a consequence of the rotational 
symmetry the light emission of microdisks is isotropic. 
For many applications, however, a directional light emis- 
sion is required. To obtain a directional emission we re- 
cently proposed to break the symmetry of a microdisk by 
placing a point scatterer inside but not at the center of 
the microdisk (see Ref. [7|). We have demonstrated that 
such a geometry leads to a significant enhancement of the 
directivity of some TM modes in outgoing light while pre- 
serving their high Q-factors. Other attempts to breaking 
the symmetry include the introduction of some other de- 
fects inside the microdisk like a linear defect [HI [3] or a 
hole [TU]. Another approach is to deform the boundary 
of the cavity [TTJ [T^ [THl [H] , or to couple light into and 
out of a microdisk with the aid of an optical fiber taper 
waveguide [15'. However, the advantage of our method 
is the analytic tractability which allows for a system- 
atic optimization of the design parameters (location and 
strength of the scatterer) with only modest numerical 
efforts. 

The purpose of this paper is to extend the theory in [7] 
to TE modes, and give a systematic study of the appear- 
ance of both highly directional TM and TE modes and its 
dependence on the distance of a point scatterer from the 
disk center. In particular, we provide arguments based 
on geometric optics to explain this dependence. This 
paper is organized as follows. In the following section. 
Sec. II, we use a Green's function method to calculate 
the positions of the resonant modes of a microdisk with 
a point scatterer in the complex wavenumber plane. In 
Sec. Ill we discuss in some detail the physical interpreta- 
tion of the parameter that describes the strength of the 
point scatterer. In Sec. IV we investigate in detail the di- 
rectivity of the resonance modes for the microdisk with 
point scatterer, and in Sec. V we show that classical ray 
optics is a useful guide to optimizing the design param- 
eters of the system. In Sec. VI we investigate the role 
of exceptional points in our system, and we finish with 
some concluding remarks in Sec. VII. 



II. THEORY OF MICRODISK CAVITIES WITH 
A POINT SCATTERER 

Let stand for in the case of TM polarization and 
for hz in the case of TE polarization. For a homogeneous 
dielectric microdisk of radius R and effective refractive 
index n in a medium of refractive index 1, Eqs. (4) and 



(6) take the form 
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inside the microdisk (r < R) and the same form with n 
replaced by 1 outside the microdisk (r > R). The reso- 
nances are obtained by imposing outgoing boundary con- 
ditions at infinity, i.e. we require that (r) cx e^^"^ 
r — > 00. Moreover, for physical reasons the value of the 
EM field at the disk center must be finite. At the bound- 
ary of the microdisk, r = R, the electric field component 
Ez and its derivative has to be continuous for TM modes. 
Similarly, for TE modes, the field hz multiplied by the re- 
fractive index and its derivative divided by the refractive 
index has to be continuous at r = i?. These boundary 
conditions lead to the 'whispering gallery' (WG) modes 
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where for TM modes, the complex wavenumbers k are 
solutions of 



A™ = J.„,{knR)H.^{kR) -nJ^{knR)H^{kR) = 0, 

(10) 

and for TE modes, the complex wavenumbers k are so- 
lutions of 

A™ = 71 Jm{,knR)H^{kR) - J^{knR)Hm{kR) = 0. 

(11) 

Here Jm and are Bessel and Hankel functions of the 
first kind respectively, m = 0,1,2,... is the azimuthal 
modal index. The constants N„i are given by 

^ ^ (H,n{kR)/J^{knR) for TM modes 
" \Hm{kR)l{nJm{nkR)) for TE modes . ^ ' 

Physically, the azimuthal modal index m = 0, 1, 2, . . . 
characterizes the field variation along the disk circumfer- 
ence, with the number of intensity hotspots being equal 
to 2m. The wavenumbers, fc, are twofold degenerate for 
TO > 0, and nondegenerate for to = 0. The radial modal 
index g = 1, 2, ... will be used to label different resonances 
with the same azimuthal modal index m. For resonances 
which are relatively close to the real axis of the the com- 
plex wavenumber plane (so called 'internal' or 'Feshbach' 
resonances), the index q gives in general the number of 
intensity spots in the radial direction. As discussed in 
detail in |16j there are exceptions to this rule for some 
TE internal resonances. 

We note that for each fixed to, there exist further so- 
lutions of Eqs. ( 10|11 1 (so called 'external' or 'shape' res- 
onances) which are typically located deeper in the lower 
half of the complex wavenumber planes compared to in- 
ternal resonances (see Refs. [IZIIIH])- The external res- 
onances are very leaky (low Q-factors) and, as a result. 
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can not be directly used for lasing. However, they are 
of theoretical interest in their own right. To properly 
distinguish between external and internal resonances one 
needs to trace the resonances as a function of the refrac- 
tive index n to the values they obtain in the limit n — > cxd. 
As discussed in detail in [16] this limiting value is real for 
internal resonances, and complex (not real) for external 
resonances. The procedure of tracing the resonances in 
the complex wavenumber plane is especially important 
for TE modes for which some of the external resonances 



lie in the same domain of the complex wavenumber plane 
as the internal resonances. For an example of this phe- 
nomenon, we refer to Sec. IV. 

In [7] we derived the Green's function for the TM 
modes of a microdisk. Following Ref. [19] it is given by a 
sum over all angular harmonics e^'^^^~^°'' multiplied by 
the corresponding radial parts. Using the same method 
it is not difficult to derive the Green's function also for 
TE waves. In fact, both functions can be written as 



G™/™ (r,ro,fc) 
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(13) 



where r< (r>) is the smaller (larger) of r and tq. The 
coefhcients are = 2 if to 7^ 0, e„j = 1 if to = 0, and 



B 



TM 



~ Jm [knR) (kR) — n J„j (knR) Jm (kR) . 



C,|/^ = Hra [knR) (kR) - n [knR) (kR) , 

= n J„, (knR) J„ (kR) - J„[ (knR) J„ (kR) , 
C™ = n H,n (knR) {kR) - H„[ (knR) {kR) . 

The resonances of the microdisk are then determined by 



the poles of the Green's function i.e. by yl™^"^^ ~ 0. 
This agrees with resonance conditions (10 1 and (111. We 



note that the resonance wavcfunctions are exponentially 
increasing as r — > cx), and hence cannot be normalized. 
However any constant multiplying the resonance wave- 
functions can be fixed by comparing the wavcfunctions 
to the residues of the Green's function. 

Using methods of self-adjoint extension theory [^I21j . 
we showed in Ref. [7] that the presence of a point scat- 
terer, which is located at a point d on the x-axis {tp = 
0), leaves the resonances of the unperturbed disk (WG 
modes) with the angular part sin(rn(f) unchanged, while 
the complex wavenumbers k^cs of the resonances with the 
angular part cos(rrup) become solutions of the transcen- 
dental equation 



1 - AGrcg(d,d,fc) ^ 0. 



(14) 



Here Greg is the regularized Green's function which is ob- 
tained by subtracting the logarithmically divergent term 
ln(fco|r — ro|)/27r from the Green's function (13l in the 
limit r, ro — > d. The two parameters A and ko can 
be absorbed in a single new parameter a defined by 
27r = — Alnfeofl- Then the condition for resonances (14) 



becomes 



= -T + :^(ln^^+7) 



27r 



4 „ ^m(^rcs) 
m— 



(15) 



where 7 0.5772 is the Euler-Mascheroni constant. The 
parameter a then determines the "strength" of the point 
scatterer located at the distance d = |d| < i? from the 
center. 

The corresponding resonance wavefunction is given by 



^'(r) =iVG(r,d,fcrcs): 



(16) 



where G is the Green's function in Eq. (13), and is a 



normalization factor. Outside of the microdisk, i.e. in 
the region r > R, the field takes the form 



N 

2 TT kyQa 



E 



, cos (mip) (^rcs^^) 



Hfn (A^res^) ■ 



(17) 

In principle, the normalization factor can be obtained 
again from the residue of the Green's function of the per- 
turbed system (see Eq. (19) in the next section). 

From Eq. ( [T5| we see that in the limit a = and 
a CX3 we recover the resonances of the unperturbed 
microdisk (without a point scatterer). For a close to 



or infinity, the approximate solutions of Eq. (15) can be 



found perturbatively. In leading order one finds 
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where Dm = {dAm/<ik)/R is given by 

D™ = Jm{knR)H^{kR) - n^J^{knR)Hm{kR), 



D. 



,TE 



Jrn{knR)Hm (kR) - J,„ {knR)Hm{kR) 
+ {n^-l) Jm{knR)H;„{kR), 



and fc is a resonance wavenuniber of the unperturbed 
disk, with aziniuthal modal index m. 



III. 



RELATION TO A FINITE SIZE 
SCATTERER 



Before we investigate the solutions of Eq. ( 15 1 and the 



corresponding resonance wave functions numerically, let 
us discuss the physical interpretation of the parameter a. 
This parameter can be related to that of a small but finite 
size scatterer as long as this scatterer can be treated in 
the s-wave approximation. The Green's function of the 
system with a point scatterer follows from self-adjoint 
extension theory [20l [21] , and is given by 



G'(r,d,fc)I?G'(d,ro,fc) 
1 -X'G«=(d,d,fc) 



G°(r,ro,fc) = G(r,ro,fc) 



(19) 

Here G'"^(d,d, fc) is the Green's function in ([13| for both 
arguments less than R and without the term —iHQ[nk\v— 
To I) /4 (i.e. without the free Green's function). The quan- 
tity V is the so-called diffraction coefficient 



V 



27r 



i7r/2 — 7 — log(nfca/2) 



(20) 



One can easily check that the resonance wavenumbers in 



( 15 1 coincide with the poles of the Green's function ( 19 1 



The form of the Green's function in ( 19 1 is identical to 
that of a system perturbed by a small s-wave scatterer 
within the Geometrical Theory of Diffraction 
The optical theorem imposes a restriction on the diffrac- 
tion coefficient given by = — 4ImI? if fc e M. It is 
a consequence of the conservation of energy. The con- 
dition for T> is equivalent to Iml?~^ = 1/4. In order to 
find a physical interpretation for the parameter a one has 



to compare the diffraction coefficient ( 20 1 with that of a 
small but finite size scatterer in the s-wave approxima- 
tion. In fact, our parameter a was chosen in such a way 



that the coefficient ( 20 ) agrees with that of a small s-wave 



scatterer of radius a with Dirichlet boundary conditions 
^ = along its circumference (see "5?]). 

Let us now perturb the microdisk by a small hole of 
radius b placed at a distance d from the center, and filled 
with a material of refractive index rif,. The diffraction 
coefficient depends only on local properties near the per- 
turbation and can be obtained from the Green's function 
for a small circular disk of radius h with refractive index 
Ub embedded within material of refractive index n (i.e. 
one can neglect the outside region with refractive index 
1). But this Green's function can be obtained in exactly 



the same manner as the Green's function in Eq. (13 1 



The only differences are the corresponding radius and 
the refractive indices. In particular, the new coefficients 
are 



A™ = nJm{knbR)Hm{knR) - 
= "^b Jm{knbR)Hm{knR) 

B™ = nJm{knbR)Jm{knR) - 
= Hb Jm{knbR)Jm{knR) - 



nb Jm{knbR)Hm{nkR), 

- n Jm{knbR)Hm{nkR), 
Hb Jm{knbR)Jm{nkR), 

- n Jm{knbR)Jm{n'kR). 



We are interested in the Green's function outside the 
small disk. If that disk is small enough to be treated 
in the s-wave approximation we need to keep only the 
term to = in the sum over modal indices 



G''(r,ro,fc) 



--Ho{kn\r 



rol) 



4 A 



Hoikn\T\)Ho{kn\ro\) . (21) 



This function has the form of ([19| with G(r,ro,fc) — 
—iHo{kn\r — ro|)/4, G^'^(r,ro,fc) = 0, and we can read 
off the diffraction coefficient as P = —HBo/Ao- One 
can check that the optical theorem is satisfied. Using 
the asymptotic forms of Bessel and Hankel functions for 
small arguments, we obtain that w i/4— [wb'^ k"^ {nl — 
n'^)]^^. The comparison with the diffraction coefficient 



for a point scatterer ( 20 1 gives the condition on the pa- 



rameter a of the point scatterer for modelling the finite 
size scatterer 



nka 
In— +7: 



(22) 



In this equation fc is real, i.e. its imaginary part which 
is small for resonances with high Q factor is neglected. 
Finally, we should note that the s-wave approximation 
itself is valid if the small but finite size scatterer is located 
not too close to the circumference of the disk, i.e. for 
\nk{R - d)| > 1, and if 



IV. FAR FIELD DIRECTIVITY 

In order to quantify the far-field behavior of the field 
^ we consider its asymptotics for r 00 which has the 
form 

^ (r, fc,e.) = ^ (r, ^, fc,es) oc ^M^f'^osr) 

Then, the directionality of the emission can be character- 
ized with the directivity D which is defined as the ratio 
of the power emitted into the main beam direction (^max 
to the total power averaged over all possible directions, 

27r |/max('^ max J I 



D 







(25) 
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Note that the resonances of the disk without a scatterer 
have D = 1 for m = and D = 2 for m 0. 

For the disk with a point scatterer, we investigate the 
level dynamics of the resonances as a function of the pa- 
rameter a by numerically solving Eq. ( 15 1 for TM and 
TE modes, respectively. As mentioned in Sec. II and 
discussed in more detail in [7] , the perturbed resonances 
reduce to the unperturbed resonances (microdisk with- 
out scatterer) in the limits a — > and a — > cxd. The 
wavenumbers of the perturbed disk thus evolve a long 
line segments in the complex wavenumber plane which 
start and end at unperturbed resonances when a is var- 
ied from to infinity. In our numerical procedure which 
consists of a Newton method to solve Eq. ([Ts]) we vary a 
from 10~'^° to 10^° and use Eq. (18 1 to obtain starting 
points for the numerics in these limits. For each res- 
onance wavenumber found this way, the corresponding 
wavefunction is obtained from Eq. (16). In particular. 



(31,1}[23.3> [20.4} (32,1) (2*,3) (TH.SJ (28.2) (33,1) (2S,3) (29,2) (22,4) 



in the outside region the wavefunction takes the form of 
(17 1 from which we can compute the directivity D. 
In fact, because of the relation between wavefunction 



and the Green's function ( 16 1 we can formally associate a 



directivity D to any point in the complex kR plane (for 
a fixed choice of d). This directivity has only a direct 
physical interpretation if the kR value corresponds to a 
resonance of the disk with a point scatterer. For gen- 
eral k, D{k) is rather a characterization of the Green's 
function (13 1. It is helpful to consider also this formal 



directivity because it enables one to clearly identify the 
regions in the kR plane where a perturbation of the cir- 
cular disk by a scatterer (at position d) should lead to 
highly directional modes. 

As a first example, we study the level dynamics of TM 
resonances in the range 12 < Re(/ci?) < 13 for a disk of 
effective refractive index n = 3.0 and radius R = 1 ^m 
(these parameters are close to the ones in Ref. [53]) with 
a point scatterer placed at two different distances (d = 
0.9 jim and d = 0.5 /um) from the center of the disk. The 
results are shown in Fig. [T] The background color in 
Fig. [l] indicates the directivity D for any point of the 
kR plane computed as mentioned above. Low values of 
D correspond to the color blue, and high values of D 
correspond to deep red. Superimposing the curves of 
the wavenumber level dynamics we can immediately see 
in which region of the kR plane the highly directional 
modes are located. 

Figure [2] shows the corresponding plot for TE modes in 
the same parameter range as in Figure [T] Like in the case 
of TM polarization there are highly directional TE modes 
for a wide range of Q-factors. From the figures one sees 
that for both TM and TE resonant modes, placing the 
point scatterer at a distance d = 0.5 fim. from the center 
yields a better directivity than for d = 0.9 /zm. We will 
explain this observation in Sec.|V] The directivity reaches 
values as high as D « 15 for some specific TM modes and 
£) « 13 for some specific TE modes. 

In Fig. |3] we show the field intensity at the distance 
r = 500 /im (far-field region) for a highly directional TM 




12.4 12.6 12,S IJ 

(31,I)<Z3,3) (20,4) (32,1) (24,3) (1I.S) (2S,2) (33,1) (1S,3) (29,2) (22,1) 

... .... . ... 



-0.04 



a: 

E 



-0.0S 



-0.12 




12.4 12.6 

Re (kR) 



2-S li 



FIG. 1: (color online). Level dynamics (black curves) of the 
TM resonances in the complex wave number plane for a di- 
electric disk with n — 3.0 and R — 1.0 /im and a point scat- 
terer of varying coupling parameter a. The solid circles mark 
the unperturbed resonances with azimuthal and radial modal 
indices (m,q). For the upper panel the scatterer is placed 
at distance d = 0.9 /zm from the centre, for the lower panel 
d = 0.5 /.im. The color code in the kR plane indicates the 
directivity D of the emission as explained in the text (blue 
marks small values of D, red marks high values of D). 



resonance mode which is obtained by a small perturba- 
tion of the unperturbed TM resonance mode with modal 
indices (21,4). It has directivity {D = 13) and complex 
wavenumber kR = 12.54929 — i 0.000045, and is com- 
pared to the unperturbed mode with directivity D = 2 
and complex wavenumber kR = 12.54876 — i 0.000001. 
The highly directional mode is obtained from the unper- 
turbed mode if a very weak point scatterer of strength 
a = 10~^ is placed at a distance d = 0.5 /xm from the 
centre of the disk. Indeed, we expect this from Fig. [l] 
because the unperturbed mode lies in a highly red region 
in the kR plane where a small perturbation should lead 



to a highly directional mode. According to Eq. (22) the 



perturbation is comparable to that of a finite size scat- 
terer of radius b w O.Olfim and refractive index ni, = I. 

In Fig. |4] we show the far field intensity of the anal- 
ogous TE mode. It is again obtained by a perturba- 
tion of the mode (21,4) of the circular disk by plac- 
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FIG. 2: (color online). Analogue of Fig. [T] for TE resonances. 
The parameters are the same as in Fig. fll 

ing a point scatterer of strength a = 10^^ at distance 
d = 0.5 /im from the centre of the disk. For this po- 
larization the unperturbed mode (21,4) is located at 
kR = 12.90089 - i 0.000001 in the complex wave num- 
ber plane while the perturbed mode is located at kR — 
12.90187 - i 0.000087. Their directivities are D = 2 and 
D = 11 respectively. Both of the above highly directional 
modes are good candidates for lasing since their Q- factors 
are large. In fact, Q « 1.4 • 10^ for the highly directional 
TM resonance and Q ~ 7.4 • 10'* for the highly direc- 
tional TE resonance. Moreover, they have unidirectional 
emission at an angular direction of 180 degrees. 

As a second example, we study the level dynamics of 
TM and TE resonances in the range 37 < Re(fci?) < 
38 for a disk of effective refractive index n = 1.4 and 
radius R = 10.51 mm (these parameters are close to the 
ones of the experimental setup of Schwefel and Preu, [25] ) 
with a point scatterer placed at two close distances d = 
9.585 mm and d = 9.385 mm from the centre of the 
disk, see figs. |5] and [6] From these two figures one can 
see that even a relatively small change in the position 
of a point scatterer can lead to a significant change in 
the level dynamics. The directivity of some perturbed 
TM modes reaches values as high as I? = 10 for a point 
scatterer located at d = 9.385 mm, see the lower panel 
of Fig. [5] and the directivity of some TE modes reaches 
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FIG. 3: The far-field intensities of the unperturbed TM res- 
onance with modal indices (21,4) (thin curve), and of the 
highly directional TM resonance (thick curve) that is ob- 
tained from the unperturbed mode by perturbing a dielectric 
disk with n = 3.0 and R — 1.0 nm with a point scatterer of 
strength a — 10"^ placed at distance d = 0.5 /im from the 
disk center. 
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FIG. 4: The field intensity in the far-field of the unperturbed 
TE resonance with modal indices (21,4) (thin curve), and of 
the highly directional TE resonance (thick curve) that is ob- 
tained from the unperturbed mode by perturbing a dielectric 
disk with n = 3.0 and R = 1.0 ^m with a point scatterer of 
strength a — 10^^ placed at distance d = 0.5 /xm from the 
disk center. 



values as high as I? = 8 for a point scatterer located at 
d = 9.585 mm, see the upper panel of Fig. [6] In general, 
however, the directivity is lower than for n = 3 in Figs, [l] 
and|2] This will be explained in the next section. 

Another interesting feature is the appearance of the 
unperturbed shape (external) TE resonance denoted by 
(23, s) in Fig. |6] which is located in the region of the 
complex wave number plane where one would expect to 
have internal resonances only. This is related to a pecu- 
liar behaviour of some of the unperturbed external TE 
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FIG. 5: (color online). Level dynamics of the TM resonances 
in the complex wave number plane for a dielectric disk with 
n = 1.4 and R — 10.51 mm and a point scatterer of varying 
coupling parameter a. The solid circles mark the unperturbed 
resonances with azimuthal and radial modal indices (m,q). 
For the upper panel the scatterer is placed at distance d = 
9.585 mm from the centre, for the lower panel d — 9.385 mm. 
As before, the color code indicates the directivity D of the 
emission. 



resonances that was observed in Ref. (TB]. For large re- 
fractive indices n, the external resonances are located 
much deeper in the complex wavenumber plane than the 
internal resonance. This clear separation by the magni- 
tude of the imaginary part of the wavenumbers ceases to 
exist for small n in the case of TE modes where some ex- 
ternal resonances mix with the internal resonances. The 
wavefunctions of such external resonances acquire similar 
features as the wavefunctions of internal resonances. In 
particular they can spoil the interpretation of the radial 
modal index q as the number of peaks of the wavefunction 
in the radial direction (see [16]) for internal resonances 
with the same azimuthal modal index. Furthermore for 
the disk with a point scatterer, the external resonance 
(23, s) in Fig. 6 illustrates the fact that unperturbed ex- 
ternal resonances also serve as starting and end points of 
the line segments that result from the level dynamics of 
perturbed resonances in the complex wavenumber plane 
upon varying a from to cxd. 

In Fig. [7] we show the field intensities at the dis- 
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FIG. 6: (color online). Analogue of Fig.jslfor TE resonances. 
The parameters are the same as in Fig.Tm The label (23, s) 
indicates a shape resonance with azimuthal model index m = 
23. 



tance r — 5000 mm (far-field region) for the unper- 
turbed resonant TM mode with modal indices (31,5), 
complex wavenumber kR — 37.599462 — i 0.488553, and 
directivity D — 2 as well as for the highly directional 
(D = 10) TM resonant mode with complex wavenumber 
kR = 37.621007 - i 0.523645. The unperturbed mode 
(31,5) transforms to the highly directional mode with 
the above complex wavenumber if we place a point scat- 
terer of strength a = 0.07 at distance d = 9.385 mm from 
the centre of the disk. We note that despite of the high 
directivity this mode is not suitable for lasing because it 
has only a small Q factor of about 36. 

In Fig. |8]we show the field intensities at the distance 
r = 5000 mm (far-field region) for the unperturbed 
resonant TE mode with modal indices (46,1), complex 
wavenumber kR = 37.129055 — 10.000177, and direc- 
tivity D — 2 'AS well as for the directional {D ^ 6) 
TE resonant mode with complex wavenumber kR — 
37.142373 - i 0.001832. The unperturbed mode (46, 1) 
transforms to the directional mode with the above com- 
plex wavenumber if we place a point scatterer of strength 
a = 4 • 10~^ at distance d — 9.585 mm from the centre 
of the disk. This mode has Q = 10* and, therefore, is 
suitable for lasing. 
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FIG. 7: The far-field intensities (r = 5000 mm) of the un- 
perturbed TM resonance with kR = 37.599462 - i 0.488553, 
m = 31, g = 5 (thin curve), and of the highly directional 
perturbed TM resonance with kR = 37.621007 - i 0.523645, 
d — 9.385 mm, a = 0.07 (thick curve) in a dielectric disk with 
n = 1.4 and R = 10.51 mm. 
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FIG. 8: The far-field intensities (r = 5000 mm) of the un- 
perturbed TE resonance with kR = 37.129055 - 10.000177, 
m = 46, q = 1 (thin curve), and of the directional per- 
turbed TE resonance with kR = 37.142373 - 10.001832, 
d = 9.585 mm, a = 4 • 10~^ (thick curve) in a dielectric disk 
with n = 1.4 and R = 10.51 mm. 

V. DIRECTIVITY AND GEOMETRIC OPTICS 

To systematically study the appearance of highly di- 
rectional modes, we calculate the average directivity for 
a region in the complex wavenumber plane as a function 
of the distance, d, of the point scatterer from the center. 
To this end we define the average directivity 
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fc,7 



fc+ 



dfc, D (fe) , (26) 



where [fc^ , k^] x [k^ , k^] is a rectangular region in the 
complex wavenumber plane of side lengths Afc^ = k^ — 



k~ and Aki = k^ — k~ . Note that the integration is over 
all k in the rectangular region where D (fc) is formally 
defined using formulas (17), (24 1, and ([25]), for general 
k rather than just for the resonant fcj-cs (see Sec. IV). In 
Figs. [9] and [TO] we show _Dav computed for the region 
12 < Re(A:i?) < 13, -0.1 < Im(fci?) < as a function 
of d for four microdisks of radius R = \ /im and effec- 
tive refractive indices of n = 3.0, n = 2.6, n = 2.25, 
and n — 1.4, and for the region 37 < Re(A:i?) < 38, 
— 1 < Im(fci?) < for four disks of radius R = 10.51mm 
with the same effective refractive indices, respectively. 
Remarkably, a rough approximation of the values for d 
which lead to high D^^ can be found from geometric op- 
tics. To this end let us consider parallel rays that come 
in from infinity and enter a dielectric disk of radius R 
and effective refractive index n. There is one ray which 
goes through the center of the disk, the central ray. The 
rays that are infinitesimally close to this central ray will 
cross the central ray at the point with distance 



R 



ri — 1 



(27) 



to the center of the disk located on the opposite side 
of the center of the disk. So, conversely, putting a point 
scatterer at this focal point leads to a strongly directional 
light emission in the an angular direction of 180° which 
also agrees with the observation in the previous section. 

The value of do is indicated by the vertical lines in 
Figs. [9] and Fig. 10 One can see that in most cases do is 
close to the optimal value range for the point scatterer 
positions in the figures. Taking into account the finite 
size of the disk, we note that formula ( 27 ) is valid only for 
refractive indices greater than 2. For smaller refractive 
indexes the optimal position should be as close as possible 
to the boundary of the disk. 



VI. EXCEPTIONAL POINTS 

The line segments that connect unperturbed reso- 
nances as the parameter a varies from zero to infinity 
can change considerably if the distance, d, of the point 
scatterer to the center is changed. This can be seen al- 
ready in Figures [5] and |6] which are for two different but 
close values of d. The connections between the unper- 
turbed resonances are very different there, i.e. the line 
segments connect different unperturbed resonances if d 
is varied only slightly. In the present section we want to 
look at this in more detail. 

As a first example we investigate the perturbation of 
one particular resonance of the microdisc. We choose 
the TM resonance with modal indices (34, 4) of the di- 
electric disk with n — 1.4 and R = 10.51. We then place 
a point scatterer at distance d from the center of the 
disk. For different values of d we vary the strength of 
the scatterer a from to infinity. This yields a family of 
line segments in the complex wavenumber plane which all 
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FIG. 9: The average directivity, Dav, of TM (thick curves) 
and TE (thin curves) polarized hght in the range 12 < 
Re{kR) < 13, —0.1 < Im(fc_R) < 0.0 for each position, d, 
of a point scatterer in microdisks of radius R = 1 /xm and 
various effective refractive indices as shown. 



FIG. 10; The average directivity, Dav, of TM (thick curves) 
and TE (thin curves) polarized light in the range 37 < 
Re(fc_R) < 38, — 1 < Im(fci?) < 0.0 for each position, d, of 
a point scatterer in microdisks of radius R = 10.51 fim and 
various effective refractive indices as shown. 



start at the unperturbed resonance (34,4) but, depend- 
ing on d, end at different unperturbed resonances. This 
is shown in Fig. [TT] It is remarkable that the connec- 
tions between the different unperturbed resonances de- 
pend so sensitively on the value of d. In the small range 
from d = 9.400 to d = 9.411 shown in Fig. [TT] the un- 
perturbed resonance (34, 4) is connected to five different 
unperturbed resonances. 



Next we investigate how the connections are rear- 
ranged. The connections between the unperturbed res- 
onances can only change upon variation of d if for one 
value of d different line segments intersect tangentially 
at a point in the complex A:i?-plane. This point corre- 
sponds to a degenerate resonance. This mechanism is 
illustrated in Fig.[l2] For d = 9.409 the resonance (34, 4) 
is connected to the resonance (9, 13), while the resonance 
(19, 9) is connected to (38, 3). For the value d = 9.411 the 
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FIG. 11: (color online) . Level dynamics of the TM resonances 
in the complex wave number plane for a dielectric disk with 
n = 1.4 and R — 10.51mm and a point scatterer of varying 
coupling parameter a for several positions d of the scatterer. 
All start from the unperturbed resonance (34,4). The solid 
circles mark the unperturbed resonances with modal indices 
(m, q). The color code indicates the directivity D of the emis- 
sion (green marks small values of D, red marks high values of 
D). 
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FIG. 12: (color online) . Level dynamics of the TM resonances 
in the complex wave number plane for a dielectric disk with 
n = 1.4 and R = 10.51mm and a point scatterer of varying 
coupling parameter a for two positions d of the scatterer. The 
solid circles mark the unperturbed resonances with azimuthal 
and radial modal indices (m,q). As the position of the scat- 
terer changes from d = 9.409 to d = 9.411 the connections of 
the unperturbed resonances change. 



connections have changed. The resonance (34, 4) is now 
connected to (38, 3), while (19, 9) is connected to (9, 13). 
Although not shown in the figure, there is a value of d 
for which two line segments touch each other at a point 
at which two perturbed resonances coallesce and become 
degenerate. 

For an open system like the dielectric disk degeneracies 
are generically of a special type which are called excep- 
tional points [23 128] . Exceptional points can be observed 
if at least two real-valued parameters of a non-hermitian 
operator are varied. In our case the non-hermitian op- 




FIG. 13: (color online). Two almost degenerate TM reso- 
nance states for a dielectric disk with n = 1.4 and R = 
10.51mm. The corresponding wavenumbers are shown by 
crosses in Fig. [12] 



erator is the differential operator acting on in Eq. ([8| 
which is non-hermitian due to the outgoing boundary 
condition. The two real parameters are the parameters 
a and d. At an exceptional point two (or more) eigen- 
values of the non-hermitian operator coallesce, and the 
corresponding eigenstates become identical. The pair of 
eigenstates which becomes degenerate at an exceptional 
point also show a characteristic behavior if the two pa- 
rameters are changed along a closed loop about the ex- 
ceptional points. Smoothly following the pair of eigen- 
states from a starting point of the loop upon one full 
traversal of the loop yields a pair of eigenstates in which 
the eigenstates started with are swapped and in which 
one member of the pair has a reversed sign. Exceptional 
points have recently received a lot of attention. For an 
example in the field of microlasers and more references 
see Ref. m\. 
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FIG. 14: Far field of the two resonance states shown by crosses 
in Fig. [12] 

To illustrate the coalescence of resonance states that 
is typical for exceptional points we show in Fig. [13] the 
fields of the two almost degenerate resonances marked 



in Fig. 12 Fig. 14 shows the corresponding far-field be- 
haviour. 

The numerical results in this section indicate that ex- 
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ceptional points are quite common for dielectric disks 
with point scatterers. They control the connections be- 
tween the unperturbed resonances and can be easily 
found by noting when these connections change upon 
varying the position d of the scatterer. 

VII. CONCLUSION 

We have shown that perturbations of a dielectric disk 
by a point scatterer can lead to highly directional reso- 
nance modes with large Q-factors. This is demonstrated 
in particular by the unidirectional modes in Figs. |3]and 
|4] To obtain modes with these properties is one of the 
main goals in the design of dielectric microcavities. 

The system studied has the advantage that it is rela- 
tively simple and can be treated to a large extent ana- 
lytically by a Green's function method. This allows for 
a systematic investigation of the system over a large pa- 
rameter range with only moderate numerical effort. We 
found that several numerical results can be understood 
with the help of a simple geometrical optics model. This 
model helps to find the optimal position of the scatterer, 
it explains why the directivity is in general higher for re- 
fractive index n = 3 than for n = 1.4, and also why the 
emission occurs predominantly in a direction opposite of 
the position of the point scatterer. It also suggests future 



investigation of an elliptical microcavity with a scatterer 
at one of the foci and refractive index inverse to the ec- 
centricity, as in that case focusing in geometric optics 
is exact, i.e. the paraxial approximation is not required 

m- 

The Green's function method also allows one to asso- 
ciate a directivity with different regions of the complex 
kR plane. This is very useful because it indicates in which 
regions of the kR plane one can expect highly directional 
modes if one perturbs the dielectric disk by the scatterer. 
It would be helpful to find a semiclassical explanation for 
the dependence of the directivity on the wavenumber k. 

Most interesting for applications we also discussed how 
the system studied can be realized physically by a small 
but finite sized scatterer. This connection can be made 
as long as the scatterer can be treated in the s-wave ap- 
proximation, and is found to be valid for examples with 
high Q-factor and directivity. An important open ques- 
tion is how the directivity and Q-factor depend on the 
size and shape of a larger scatterer when corrections to 
the s-wave approximation are taken into account. 
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